Orbital Degree of Freedom and Phase Separation 
in Ferromagnetic Manganites at Finite Temperatures 
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The spin and orbital phase diagram for perovskite manganites are investigated as a function of 
temperature and hole concentration. The superexchange and double exchange interactions dominate 
the ferromagnetic phases in the low and high concentration regions of doped holes, respectively. 
The two interactions favor different orbital states each other. Between the phases, two interactions 
compete with each other and the phase separation appears in the wide range of temperature and 
hole concentration. The anisotropy of the orbital space causes discontinuous changes of the orbital 
state and promotes the phase separation. The relation between the phase separation and the stripe- 
and sheet-type charge segregation is discussed. 



PACS numbers; 75.30.Vn, 75.30.Et, 71.10.-w 
I. INTRODUCTION 

Doped perovskite manganites and their related com- 
pounds have attracted much attention, since theyj^how 
not only the colossal magnetoresistance (CMRp Q but 
many interesting phenomena such as a wide variety of 
magnetic structure, charge ordering and structural phase 
transition. Although the ferromagnetic phase commonly 
appears in the manganites, the origin still remains to be 
clarified. Almost a half-century ago, the double exchange 
(DE) interaction was proposed to explain the close con- 
nection between the appearance of ferromagnetism and 
the metallic conductivity.^ In the scenario, the Hund 
coupling between carriers and localized spins is stressed. 
It has been recognized that the ferromagnetic metal- 
lic state in the highly doped region of Lai-jjA^jMnOa 
{x ~ 0.3) with A being a divalent ion is understood based 
on this scenano, where the compounds show the wide 
band width. QQ 

On the contrary, thej-DE scenario is not applied to 
the lightly doped regionlj {x < 0.2) where the CMR ef- 
fect is observed. In the region, the degeneracy of eg 
orbitals in a Mn'^^ ion exists and affects the physical 
properties. The degeneracy is called the orbital degree 
of freedom. With taking into account the orbital de- 
gree together with the electron correlation, the additional 
ferromagnetic interaction, that is, the ferromagnetic su- 
perexchange (SE) interaction, is derived. This is associ- 
ated with the alternate alignment oLlha orbital termed 
antiferro(AF)-type orbital ordering.Ehy The SE inter- 
action dominates the ferromagnetic spin alignment ob- 
served in the a6-plane in LaMnOa and the quasi fwp-i 
dimensional dispersion relation of the spin wave in it.E3Ell 
When holes are introduced into the insulating LaMnOa, 
the successive transitions occur in magnetic and trans- 
port phase diagrams; with increasing x, it is observed in 
Lai_a;Sra;Mn03 as almost two dimensional ferromagnetic 
(A-type AF) insulator — > isotijopic ferromagnetic insu- 
lator — > ferromagnetic metal.clEj The first order phase 
transition between two ferromagnetic states recently dis- 



covered in Lai_j;Sra;Mn03 with x ~ 0.1313 indicates that 
the orbital state also changes at the transition. In order 
to understand a dramatic change of electronic states in 
lightly doped region and its relation to CMR, it is indis- 
pensable to study the mutual relation between the two 
ferromagnetic interactions, i.e., DE and SE. 

In this paper, we investigate the spin and orbital phase 
diagram as a function of temperature (T) and hole con- 
centration (x). We focus on the competition and cooper- 
ation between the two ferromagnetic interactions SE and 
DE. We show that the SE and DE interactions dominate 
the ferromagnetic phases in the low and high concentra- 
tion regions of doped holes, respectively, and favor the 
different orbital structures each other. Between the two 
phases, the phase separation (PS) appears in the wide 
range of x and T. It is shown that the phase separation 
is promoted by the anisotropy in the orbital space. The 
spin and orbital phase diagram at T = was obtained by 
the Hartree-Fock theory and interpreted in terms of the 
SE and DE interactions in Ref. |l^. The PS state between 
two ferromagnetic phases driven by the DE interaction 
and the Jahn- Teller distortion at T = was discussed in 
Ref. |lj. In this paper, we obtain the PS state based on 
the model with strong correlation of electrons at finite T. 

In Sect. II, the model Hamiltonian, where the electron 
correlation and the orbital degeneracy are taken into ac- 
count, is introduced. In Sect. Ill, formulation to calculate 
the phase diagram at finite T and x is presented. Numer- 
ical results are shown in Sect. IV and the last section is 
devoted to summary and discussion. 



II. MODEL 

Let us consider the model Hamiltonian which describes 
the electronic structure in perovskite manganites. We set 
up the cubic lattice consisting of manganeses ions. Two 
Cg orbitals are introduced in each ion and t2g electrons are 
treated as a localized spin {St^ ) with S — 3/2. Between 
e„ electrons, three kinds of the Coulomb interaction, that 



is, the intra-orbital Coulomb interaction (U), the inter- 
orbital one {U') and the exchange interaction(/), are 
taken into account. There also exist the Hund coupling 
(Jh) between eg and t2g spins and the electron trans- 
fer tj^ between site i with orbital 7 and site j with 7'. 
Among these energies, the Coulomb interactions are the 
largest one. Therefore, by excluding the doubly occupied 
state at each site, we derive the effective Hamiltpaian de- 
scribing the low energy spin and orbital states :0 



Ti. — Tit + Ti-j + Ti-H + Haf ■ 



(1) 



The first and second terms correspond to the so-called 
t- and J-terms in the tJ-model for Cg electrons, respec- 
tively. These are given by 



(ij)77'cr 



■J7'cr 



+ H.C 



(2) 



manganites, as well as the undoped insulator. 2) Since 
Ji > J2 , the ferromagnetic state associated with the AF- 
type orbital order is stabilized by Tij. Therefore, two 
kinds of the ferromagnetic interaction, that is, SE and 
DE arc included in the model. 3) As seen in Hj, the or- 
bital pseudo-spin space is strongly anisotropic unlike the 
spin space. 



III. FORMULATION 

In order to calculate the spin and orbital states at 
finite temperatures and investigate the phase separa- 
tion, w e g eneralize the mean field theory proposed by de 
GennesO Hereafter, the spin (5) and pseudo-spin (T) 
variables are denoted by u in the unified fashion. In this 
theory, the spin and orbital pseudo-spin are treated as 
classical vectors as follows 
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where 



/27r , 
cos! — n; ]Ti, 



/2n . 
smi — n; )T„ , 
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and (rix , riy 



di~f(j is the annihila- 



and 



^H — —Jh 2_^ ^t2gi ■ ^i J 



HaF — JaF 2_^ ^t2gi • St2g] 



(5) 



(6) 



The detailed derivation of the Hamiltonian is presented 
in Ref. O. Main features of the Hamiltonian are sum- 
marized as follows: I) This is applicable to the doped 
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9f cos (j)f , sin et sin (j)f , cos 6','' ) , (7) 



and 



(3) {T-,Ty,Tn^^{sm0l,O,cos9l) 



(8) 



(1,2,3). / denotes the direction of 



bond connecting i and j sites. 
tion operator of Cg electron at site i with spin a and 
orbital 7 with excluding double occupancy. St is the 
spin operator of the eg electron and Ti is the pseudo- 
spin operator for the orbital degree of freedom defined as 
T^ = (1/2) E.77' ^7-('^)77'rfnV. Ji = WiU' - I) and 
J2 = t§ /([/' + / + 2 J If) where io is the transfer intensity 
between d^z^_j.2 orbitals in the z-direction, and the rela- 
tion U = U' + I \s assumed. The orbital dependence of 

t]7 is estimated from the Slater-Koster formulas. The 

__ 

third and fourth terms in Eq. (pi) describe the Hund cou- 
pling between Cg and t2g spins and the antiferromagnetic 
interaction between t2g spins, respectively, as expressed 
as 



where the motion of the pseudo-spin is assumed to be 
confined in the xz-plane. 6\ in Eq. (ph characterizes the 
orbital state at site i as 

1^*) = cos{e\/2)\d^,2_,2) + sin[e\/2)\d^2^y2) . (9) 

t2g spins are assumed to be parallel to the Cg one. The 
thermal distributions of the spin and pseudo-spin are de- 
scribed by the distribution function which is a function 
of the relative angle between Ui and the mean field A" , 



K{u,) = —e^]){Xi ■m'i) , 



(10) 



where 7n^{= Ui/\u\) is termed the spin (pseudo-spin) 
magnetization and the normalization factor is defined by 



I dO j d(/)exp(A"cos6') , 



and 



2ir 



(i0exp(A'cos6'). 



(11) 



(12) 



The mean fields are assumed to be written as A" — 
A"(sin0", O,cos0^). By utilizing the distribution func- 
tions defined in Eq. (lOh, the expectation values of oper- 
ators Ai{S) and Bi{T) are obtained as 



(A.) 
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and 



{B. 



2tt 



de'wl{fi)B(f) , 



(14) 



respectively. In this scheme, the free energy is repre- 
sented by summation of the expectation values of the 
Hamiltonian and the entropy of spin and pseudo-spin as 
foUows: 



T^{n)-NT{S' +S') . 



(15) 



N is the number of Mn ions and 5" is the entropy cal- 
culated by 



5" = -(lnu;"(u))„ 



(16) 



By minimizing J-' with respect to A" and Qi, the mean 
field solutions are obtained. It is briefly noticed that the 
above formulation gives the unphysical states at very low 
temperatures (T < Tneg ~ Ji{2)/^0) where the entropy 
becomes negative. Therefore, we restrict our calculation 
in the region above Tneg- However, at T = 0, the spin 
and orbital states are calculated without any trouble in 
the entropy with the assumption of the full polarizations 
of spin and pseudo-spin. 

Next, we concentrate on the calculation of (H) in 
Eq. ( p^ . As shown in Eq. (^) , Tij is represented by S and 
T. By introducing the rotating frame in the spin(pseudo- 
spin) space, the z-component of the spin (pseudo-spin) 
in the frame is given by 



= cos ei" 



sine"? 



(17) 



sin 9f sin 0^ as expected from the double exchange 
interaction.cl By diagonalizing the energy in the momen- 
tum space. Tit is given by 



Ni 

1 1=1 



(20) 



where I indicates the band of ft-^g and Ni is the number 
of the bands, e^ corresponds to the energy of the l-th 
band. As a result, the expectation value of Ht per site is 
obtained by 



Et 



t 1=1 



(21) 



which is a function of the spin and pseudo-spin angles 
at each site, {&f} and {6*}, and the amplitudes of the 
mean fields. A* and A*, ep in Eq. (|l]) is the fermi energy 
of h^^ determined in the equation, 



1 ^' 

TV ^ ^ • 
k '=1 



jFie'j^-SF), (22) 

where /f(£) is the fermi distribution function. 
IV. NUMERICAL RESULTS 



which is parallel to the mean field A". Thus, {u^)u is 
adopted as the order parameter which has the relation, 
{ui)u = 5(™"^)n- The spin part in Hj is rewritten by 
using {mf^)s and the relative angle in the spin space as 
{fh^^')1 cos(8f — 0p. On the other hand, the orbital part 
includes the term cos(8* -|-0*), which originates from the 
anisotropy in the orbital space. TLa^ is also rewritten 
by using {mf^)s and 8" under the relation of {S)s = 
4:{St2g)s- As faiLthe transfer term Tit, we introduce the 
rotating frameEl and decompose the electron operator 
as di-ya = h\zf^zl^ where h\ is a spin- less and orbital- 
less fermion operator and zf^ and z*^ are the elements of 
the unitary matrix (C/'*^")) in the spin and pseudo-spin 
frames, respectively. These are defined by 



[/" 



H] 






(18) 



i<l>t/2 



with zf^ = cos(6'f/2)e-*"^'/2 ^nd z|^ = sin(6l|/2) e*' 
for spin, and z*| = cos(0*/2) and zl^ = sin(6'*/2) for 
orbital. By using the form. Tit is rewritten as 



Ut 



E^ 



r h h^ 



H.c. 



(19) 



and i* 



with tlj = Ea^l<T*^ja: — -^J 

The former gives e'("*?-"*.°)/2 cos 6 



E..-' z^*^h' z: 



'cos 6*1 -e-^'^'I'^'-'f''^^/^ 



A. phase diagram at T = 

In this subsection, we show the numerical results at 
T = 0. For examining both spin and orbital orderings, 
two kinds of sublattice are introduced. We assume ferro- 
magnetic (F)-type and three kinds of antiferro (AF)-type 
spin (pseudo-spin) orderings, which are layer (A)-type, 
rod (C)-type and NaCl (G)-type. 

In Fig. 1(a), the ground state energy (Egs) is shown 
as a function of hole concentration (x) for several values 
of Jaf /to- Double- or multi-minima appear in the Eqs- 
X curve depending on the value of JAp/to- Therefore, 
the homogenous phase is not stable against the phase 
separation. This feature is remarkable in the region of 
0.1 < a; < 0.4. In Fig. 1(b), Egs is decomposed into 
(Ht) and (Hj) for JAp/to = 0. By drawing a tangent 
line in the Eqs-x curve as shown in Fig. 1(a), the phase 
separation is obtained. By using the so-called Maxwell 
construction, the phase diagram at T = is obtained 
in the plane of Jaf and x (Fig. 2). The parameter val- 
ues are chosen to be Ji/to — 0-25 and J2/to ~ 0.0625. 
JAp/ta for manganites is estimated from the Neel tem- 
perature in CaMnOs to be 0.001 ^ 0.01. Let us consider 
the case of JAp/to — 0.004. With doping of holes, the 
magnetic structure is changed as A-AF — > PS(A-AF/Fi) 
^ Fi ^ PS(Fi/F2) -> F2, where PS(A/B) implies the 
phase separation between A and B phases. The canted 




FIG. 1. The ground state energy [Eos) as a functic 
of hole concentration (a;), (a): JAp/to is chosen to be i 
0.004, and 0.01. The broken hnes and the filled triangl( 
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FIG. 3. The densities of state (DOS) for the spin-less and 
orbital-less fermions /ij (a) in the Fi phase and (b) in the F2 
phase. The shaded areas show the occupied state of /ig 
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FIG. 2. The phase diagram at T = in the plane of anti- 
ferromagnetic interaction Jaf and hole concentration x. Fi 
and F2 are the ferromagnetic phases with different types of 
orbital ordering. PS(Fi/F2) is the phase separated state be- 
tween the Fi and F2 phases. Types of orbital ordering in the 
two phases are schematically presented. In the dotted region, 
there exist PS(A-AF/Fi) and PS(A-AF/C-AF). The param- 
eter values are chosen to be Ji/to ~ 0.25 and J2/fo = 0.0625. 



spin structure does not appear. Fi and F2 are the two 
kinds of ferromagnetic phase discussed below in more de- 
tail. Between Fi and F2 phases, the PS state appears and 
dominates the large region of the phase diagram. For ex- 
ample, at a; = 0.2, the Fi and F2 phases coexist with the 
different volume fractions of 60% and 40%, respectively. 
We also find the PS state between A-AF and Fi phases 
in the region of 0.0 < a; < 0.03. 

Now we focus on two kinds of ferromagnetic phase and 
the PS state between them. The Fi and F2 phases orig- 
inate from the SE interaction between Cg orbitals and 
the DE one, respectively. The interactions have different 
types of ocbital ordering as shown in Fig. 2. These are 
the C-typeEl with {6^/9%) = (7r/2, 3tt/2) and the A-type 
with {6\/9g) = (7r/6, — 7r/6), respectively, where 0*^^,^-. is 
the angle in the orbital space in the A{B) sublattice. It 
is known that the AF-type orbital ordering obtained in 
the Fi phase is favorable to the ferromagnetic SE in- 
teraction through the coupling between spin and orbital 
degrees in Hj. On the other hand, the F-type orbital 
ordering promotes the DE interaction by increasing the 
gain of the kinetic energy. To show the relation between 
the orbital ordering and the kinetic energy, we present 
the density of state (DOS) of the spin-less and orbital- 
less fermions in the Fi and F2 phases in Fig. 3(a) and 
(b) , respectively. It is clearly shown that the band width 
in the F2 phase is larger than that in the Fi phase. In 
addition, DOS in the F2 phase has a broad peak around 
—2 < uj/to < —0.8 which results from the quasi-one di- 
mensional orbital ordering. Because of the structure in 
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FIG. 4. A sequential change of orbital states as a tunc 
of hole concentration x. 0\rg-. is the angle in the orbital s] 
in the A{B) orbital sublattice. The schematic orbital s 
are shown. In the phase-/ and -II, the dotted areas show 
region where the hole concentration is rich. 



DOS, the kinetic energy further decreases in the F2 pi 
more than the Fi phase. 

In order to investigate the stability of the PS state 
pearing between the Fi and F2 phases, the ground s 
energy is decomposed into the contributions from the 
interaction {{Hj)) and the DE one {{Ht)) (see Fig. 1( ,, 
We find that with increasing x, {Hj) increases and (Ht) 
decreases. Several kinks appear in the {Hj)-x and {T-it)-x 
curves, which imply the discontinuous change of the state 
with changing x. The PS(Fi/F2) state shown in Fig. 2 
corresponds to the region, where the two ferromagnetic 
interactions compete with each other and the discontinu- 
ous changes appear in the {Ti.ju\)-x curve. In Fig. 4, we 
present the x dependence of the orbital state assuming 
the homogeneous phase. It is clearly shown that the dis- 
continuous change of {'Hj(t) )-x curve is ascribed to that of 
the orbital state. In particular, in the phase-I and -II, the 
symmetry of the orbital is lower than that in the Fi and 
F2 phases and the stripe-type (quasi one dimensional) 
and sheet-type (two dimensional) charge disproportion 
is realized, respectively. These remarkable features origi- 
nate from the anisotropy in the orbital pseudo-spin space. 
We also note that because of the anisotropy, the orbital 
state dose not change continuously from Fi to F2. It is 
summarized that the main origin of the PS state in the 
ferromagnetic state is 1) the existence of two kinds of fer- 
romagnetic interaction which favor the different types of 
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FIG. 5. The phase diagram in the plane of temperature 
(T) and hole concentration (x). The homogeneous state is as- 
sumed. The straight and dotted lines show the ferromagnetic 
Curie temperature (Tc) and the orbital ordered temperature 
i!Too), respectively. The parameter values are chosen to be 
Ji/io = 0.25, J2A0 = 0.0625 and JafIU = 0.004. 

orbital state, and 2) the discontinuous change of orbital 
state due to the anisotropy in the orbital space unlike the 
spin case. 



B. phase diagram at finite T 

In this subsection, we show the numerical results at 
finite T and discuss how the PS state changes with T . 
As the order parameter of spin, we assume the ferromag- 
netic ordering and focus on the Fi and F2 phases and the 
PS state between them. We consider the G- and F-type 
orbital orderings which are enough to discuss the orbital 
state in the ferromagnetic state of the present interest. 

In Fig. 5, the phase diagram is presented at finite T 
where the homogeneous phase is assumed. Parameter 
values are chosen to be JAp/to = 0.004, Ji/io = 0.25 
and J2/io = 0.0625. At x = 0.0, the orbital ordered tem- 
perature {Too) is higher than the ferromagnetic Curie 
temperature (T^), because the interaction between or- 
bitals (3J1/2) in the paramagnetic state is larger than 
that between spins (Ji/2) in the orbital disordered state, 
as seen in the first term in Tij. With increasing x, Tq 
monotonically increases. On the other hand. Too de- 
creases and becomes its minimum around x ~ 0.25. This 
is the consequence of the change of orbital ordering from 
G-type to F-type. The G- and F-type orbital orderings 
are favorable to the SE and DE interactions, respectively, 
so that the orderings occur in the lower and higher x re- 
gions. In Fig. 6(a), we present the free energy as a func- 
tion of X at several temperatures. For T/io < 0.025, the 
double minima around x = 0.1 and 0.4 exist as discussed 
in the previous subsection at T = 0. With increasing 
T, the double minima are gradually smeared out and a 
new local minimum appears around x = 0.3. It implies 
that another phase becomes stable around x = 0.3 and 




FIG. 6. The free energy as a function of hole concentrat 
(x). (a): T/to is chosen to be 0, 0.04, and 0.15. The bro 
lines and the filled triangles show the tangent lines of the J 
curve and the points of contact between the two. (b): J 
decomposed into the contributions from TS, {Tit) and (7i 
T/to is chosen to be 0.04. 



two different kinds of the PS state appear at the tt 
peraturc. With further increasing temperature, seve 
shallow minima appear in the J--x curve. Finally, 
fine structure disappears and the homogeneous phase 
comes stable in the whole region of a;. In Fig. 6(b), 
free energy is decomposed into the contributions from 
TS, {Ht) and (Tij) at T/to = 0.04. 

By applying the Maxwell construction to the free en- 
ergy presented in Fig. 6(a), the PS states are obtained 
and presented in Fig. 7. The PS states dominate the 
large area in the x-T plane. A variety of the PS states 
appears with several types of spin and orbital states. 
Each PS state is represented by the combination of spin 
and orbital states, such as PS(spin-P, orbital-G/spin- 
F, orbital-P) for PS-III and PS(spin-F, orbital-G/spin- 
F, orbital-F)=PS(Fi/F2) for PS- VII. Here, P indicates 
the paramagnetic (orbital) state. It is mentioned that 
the phase diagram in Fig. 7 has much analogy with 
that in eutectic alloys. For example, let us focus on 
the region below T/to — 0.05. Here, the Fi and F2 
phases and PS- VII correspond to the two kinds of ho- 
mogeneous solid phases, termed A and B, and the PS 
state between them (PS(A/B)) in binary alloys, respec- 
tively. In the case of the binary alloys, the liquid(L)- 
phase becomes stable due to the entropy at high temper- 
atures. Thus, with increasing temperature, the succes- 
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FIG. 7. The phase diagram at finite temperatures. The 
shaded area shows the phase separated region. The spin 
and orbital states in each state is PS-I: PS(spin-P, orbital-G 
/spin-P, orbital-P), PS-II: PS(spin-P, orbital-P /spin-F, or- 
bital-P), PS-III: PS(spin-P, orbital-G /spin-F, orbital-P), 
PS-IV: PS(spin-P, orbital-G /spin-F, orbital-G), PS-V: 
PS(spin-F, orbital-G /spin-F, orbital-P), PS- VI: PS(spin-F, 
orbital-P /spin-F, orbital-F), and PS- VII: PS(spin-F, or- 
bital-G /spin-F, orbital-F)=(Fi/F2). The parameter values 
are the same as those in Fig. 5. 



sive transition occurs as PS(A/B) -^ (PS(L/A(B))) -^ 
L. The states, L, PS(L/A) and PS(L/B), correspond to 
the (spin-F, orbital-P) phase, PS-V, and PS- VI in Fig. 7, 
respectively. By the analogy between two systems, the 
point at T/to — 0.025 and x = 0.27 corresponds to 
the eutectic point. In the J--x curve shown in Fig. 6, 
above three states reflect on the three minima observed 
at T/to = 0.004. By decomposing the free energy into 
the three terms: {Tij), {Ht) and TS, we confirm that 
the middle part corresponding to the (spin-F, orbital-P) 
phase is stabilized by the entropy. 

In Fig. 8, we present effects of the magnetic field (B) 
on the phase diagram. The magnitude of the applied 
magnetic field is chosen to be gfisB/to — 0.02 which 
corresponds to 50 Tesla for to = O.SeV and g — 2. We 
find that the PS state shrinks in the magnetic field. The 
remarkable change is observed in PS-II and III where the 
spin-F and -P phases coexist. The magnetic field stabi- 
lizes the ferromagnetic phase so that the PS states are 
replaced by PS-V and the uniform ferromagnetic state. 
The region of PS- VII (PS(Fi/F2)) is also suppressed in 
the magnetic field. Because the magnitude of the mag- 
netization in the phase Fi is smaller than that in the 
F2 phase, the magnetic field increases the magnetization 
and stabilizes the Fi phase. 



V. SUMMARY AND DISCUSSION 

In this paper, we study the spin and orbital phase dia- 
gram for perovskite manganites at finite T and x. In par- 
ticular, we pay our attention to two kinds of ferromag- 
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FIG. 8. The phase diagram at finite temperatures in the 
applied magnetic field (B). The open and filled circles show 
the boundary of the phase separated region in qjibB /ta — 
and 0.02, respectively. The other parameter values are the 
same as those in Fig. 5. 



netic phase appearing at different hole concentrations. 
The SE and DE interactions dominate the ferromagnetic 
phases in the lower and higher x and favor the AF- and F- 
type orbital orderings, respectively. Between the phases, 
the two interactions compete with each other and the 
phases are unstable against the phase separation. The 
PS states at finite T have much analogy with that in the 
binary alloys. 

It is worth to compare PS(Fi/F2) with PS(AF/F). 
As shown in Fig. 2 (Jaf/^o = 0.004), PS(Fi/F2) ap- 
pears in the region of higher x than PS(A-AF/F). This 
originates from the following sequential change of the 
state with doping of holes as I:(spin-A,orbital-G) -^ 
II:(spin-F,orbital-G) -^ III:(spin-F, orbital-F). The or- 
bital state changes at higher x than the spin state. As 
a resuh, PS(A-AF/F) and PS(Fi/F2) appear between I 
and II, and II and III, respectively. This is because 1) 
at X = 0, the ferromagnetic interaction between spins 
is weaker than the AF one between orbitals, as men- 
tioned in Sect. IV B, and 2) at a; = 0, the AF interac- 
tion along the c-axis is much weaker than the ferromag- 
netic one in the a6-plane. We also notice in Fig. 2 that 
PS(Fi/F2) dominates a larger region in the phase dia- 
gram than PS(A-AF/F). This mainly results from the 
anisotropy in the orbital pseudo-spin space. As shown 
in Fig. 4, 0*^,^s indicating the orbital state discontinu- 
ously changes with x in the region of 0.06 < x < 0.41. 
Continuous change from Fi to F2 is prevented by the 
anisotropy in the orbital space. This is highly in contrast 
to the spin case where the incommensurate and/or flux 
states associated with the continuous change of the spin 
angle become more stable than some PS states.cj The 
anisotropy in the orbital space also stabilizes the homo- 
geneous state in the region of a; < 0.06. On the other 
hand, PS(A-AF/F) appears by doping of infinitesimal 
holes. Furthermore, the microscopic charge segregation 
appearing in the phase-I and -II (Fig. 4) is also due to the 



orbital degree of freedom. Here, the stripe- or sheet-type 
charge disproportion is realized and the SE and DE inter- 
actions dominate different microscopic regions (bonds). 
These unique phases are ascribed to the dimensionality 
control of charge carriers through the orbital orderings. 
It is mentioned that when the ocbital degree of freedom is 
taken into account, PS(AF/F)lj discussed in the double 
exchange modelEfl is suppressed. This is because A-AF 
is realized at a; = instead of G-AF and the ratio of 
the band width between A-AF and F is Waf/Wf=2/3. 
This ratio is much larger than that between G-AF and 
F which is of the order of 0{ta/JH). Therefore, the 
PS region, where the compressibility (k = {dii/dx)~^) 
is negative, shrinks. The {d^j.2_^2 / d'iy2_j.2)-iy^e orbital 
ordering expected from the lattice distortion in LaMnOa 
further enhances Waf/Wf, because the transfer inten- 
sity along the c-axis is reduced in the ordering. 

It should be noticed that the following effects may sup- 
press the phase separation discussed in the paper. In the 
present calculation, the order parameters for spin and 
orbital are restricted in a dice consisting of 2 x 2 x 2 
Mn ions. Other types of the ordering become candidates 
for the solution with the lower energy, especially, in the 
lightly doped region. However, the orbital ordering with 
the long periodicity iSjiess important in comparison with 
that in the spin case.ES The orbital ordering associated 
with continuous change of the pseudo-spin is prohibited 
by the anisotropy, as discussed above. Neither the quan- 
tum fluctuation neglected in the mean field theory nor 
the long range Coulomb interaction favor the phase sep- 
aration. When the effects are taken into account, the area 
of PS in the x-T plane shrinks and certain regions will 
be replaced by the homogeneous phases. In this case, it 
is expected that the phases with the microscopic charge 
segregation, such as the phase-I and -II shown in Fig. 4, 
remain, instead of the macroscopic phase separation. 

For observation of the PS(Fi/F2) state proposed in 
this paper, the most direct probe is the resonant x-ray 
scattering which has recently been d yrdo ped as a tech- 
nique to observe the orbital ordering.E3E3 Here, the de- 
tailed measurement at several orbital reflection points 
are required to confirm PS where different orbital order- 
ings coexist. Observation of the inhomogeneous lattice 
distortion is also considered as one of the evidence of 
PS(Fi/F2), although this is an indirect one. Several ex- 
perimental results have reported an inhomogenenity in 
the lattice degree of freedom. In Lai_a;Sra;Mn03, two 
kinds of Mn-0 bond with different lengthsp-are observed 
by the pair distribution function analyses E3 These val- 
ues are almost independent of x, although the averaged 
orthohombicity decreases with x. Since two kinds of the 
bond are observed far below Tc where the magnetiza- 
tion is almost saturated, PS(AF/F) is excluded and PS 
with different orbital orderings explain the experimental 
results. The more direct evidence of PS was reported bji 
the synchrotron x-ray diffraction in Lao.88Sro.12MnO3.Lj 
Below 350K, some of the diffraction peaks split and the 
minor phase with 20% volume fraction appears. This 



phase shows the larger orthohombic distortion than the 
major one in the region of 105K < T < 350K. Thus, the 
experimental data are consistent with the existence of 
PS(Fi/F2) where the major and minor phases correspond 
to the F2 and Fi phases, respectively. In this compound, 
the first order phase transition from ferromagnetic ip^ 
sulator to ferromagnetic metal occurs at T = 145K.t3 
Through the systematic experiments, it has been revealed 
that this magnetic transition is ascribed to the transi- 
tion between the orbital ordering and disordering. The 
experimental results strongly suggest that the two differ- 
ent interactions, i.e., SE and DE, are concerned in the 
transition and unconventional experimental results are 
understood in terms of the interactions. It is desired to 
carry out further experimental and theoretical investiga- 
tions to clarify roles of the PS state on the unconventional 
phenomena. 
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